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motion  result  from  the  roller/race  contact  resonance. 
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ABSTRACT 
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An  analytical  formulation  for  the  roller  motion  in  a cylindrical  roller 
bearing  is  presented  in  terms  of  the  classical  differential  equations  of 
motion.  Roller-race  Interaction  is  analyzed  in  detail  and  the  resulting 
force  and  moment  vectors  are  determined.  Formulation  for  the  roller  end 
and  race  flange  interaction  during  skewing  of  the  roller  is  also  considered. 
Elastohydrodynamic  traction  models  for  the  turbo-33  type  oil,  5P4E  polyphenyl 
ether,  and  the  MIL-L-7808  oil  are  sunmarized  for  the  purpose  of  computing 
roller/race  traction.  The  natdre  of  the  general  motion  of  the  roller  is 
Investigated  by  integrating  the  equations  of  motion  for  a typical  roller 
bearing.  It  is  found  that  very  high  frequency  vibrations  in  the  roller 
motion  result  from  the  roller/race  contact  resonance. 
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SECTION  I 
INTRODUCTION 

Practical  rotor  ayatams  subjected  to  relatively  large  radial  loads  are  often 
supported  by  roller  bearings  primarily  due  to  the  large  load  support  end  stiff- 
ness characteristics  of  such  bearing.  A dominant  radial  load,  however,  causes 
a number  of  rollers  to  be  completely  unloaded  as  they  travel  In  orbit.  Thus, 
each  roller  la  subjected  to  large  load  variations  and,  therefore,  the  resulting 
motion  Is  somewhat  complicated.  Other  factors  such  as  Irregular  traction  charac- 
teristics, roller-cage  Interactions,  relative  race  misalignment,  etc.,  further 
complicate  the  roller  motion  and  sometimes  lead  to  Instabilities  termed  as  "skid- 
ding" and  "skewing".  Skidding  normally  represents  large  sliding  between  the 
roller  and  the  race,  while  skewing  denotes  rotation  of  the  roller  about  Its  trans- 
verse axes.  The  primary  objective  of  a bearing  system  design  Is  to  reduce  or 
eliminate  such  Instabllltlea  under  the  prescribed  operating  conditions.  A realistic 
_ design  tool  must  simulate  the  roller  motion  with  adequate  generality  In  order  to 

B provide  the  required  design  guidelines.  The  development  of  such  an  analytical  tool 

Is  the  primary  objective  of  Che  present  investigation. 


Most  of  the  analytical  simulations  of  roller  bearing  performance  available  to  date 
have  been  restricted  to  quasi-static  type  of  force  balance  models.  It  la  clear 
chat  although  such  models  will  give  adequate  Information  about  the  overall  load 
distribution  and  bearing  fatigue  life,  they  fall  to  provide  any  Insight  Into  the 
general  motion  of  the  roller  and  the  resulting  dynamic  behavior  of  Che  bearing. 

For  Che  simulation  of  even  some  of  Che  simplest  forms  of  roller  motions.  It  Is 
necessary  chat  the  general  differential  equations  of  motion  be  formulated  and  In- 
tegrated In  time,  with  prescribed  Initial  conditions.  The  external  forces  on  the 
roller  are  determined  by  the  applied  loads  on  the  bearing,  roller-race  Interaction 
models,  lubricant  traction  behavior,  roller-cage  Interactions  and  some  ocher 
factors  such  as  lubricant  drag  and  churning. 

This  report  contains  the  analysis  relevant  to  the  above  general  Interactions  in  the 
case  of  a cylindrical  roller  bearing.  The  work  Is  somewhat  similar  Co  an  earlier 
Investigation  (reference  2)  for  the  case  of  ball  bearings.  The  motion  of  the  cage, 
which  Is  a crucial  element  for  both  ball  and  roller  bearings,  has  also  been 
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formulated  in  an  earlier  report.  Thus,  this  report  Is  the  third  part  of  a con- 
tinuing project  with  the  ultimate  objective  of  developing  generalized  computer 
simulation  models  for  the  dynamic  performance  of  both  ball  and  roller  bearings 
The  next  forthcoming  part,  IV,  will  contain  details  about  the  traction  models, 
churning  losses  and  related  factors.  Some  results  of  parametric  studies  will  be 
described  In  Part  V and  Parts  VI  and  VII  will  respectively  be  the  user's  manual 
and  compilation  listing  of  the  computer  program. 

The  generalized  equations  of  motion  for  a cylindrical  roller  are  described  In 
the  next  section  of  this  report.  Sections  III  and  IV  are  devoted  to  roller-race 
interactions  and  the  race  flange  and  roller  end  Interaction.  Simply  for  the  pur- 
pose of  computing  the  initial  conditions  for  the  general  equations  of  motions, 

Che  static  and  quasi-static  equilibrium  equations  are  contained  in  Sections  V 
and  VI,  respectively. 


SECTION  II 


GENERALIZED  EQUATIONS  OF  ROLLER  MOTION 


Similar  to  a ball  bearing,  the  motion  of  a roller  In  a generalized  simulation  of 
the  dynamic  performance  of  a roller  bearing  can  be  considered  In  two  parts: 


(1)  Motion  of  the  roller  mass  center  In  an  Inertial  reference  frame. 


(2)  Rotational  motion  of  the  roller  about  Its  mass  center. 


A formulation  of  both  the  above  parts  of  roller  motion  and  a definition  of  the 
relevant  coordinate  frames  will  be  the  basic  objectives  of  this  section. 


Motion  of  Ball  Mass  Center 


As  shown  In  Figure  II-l,  the  translational  motion  of  the  roller  center  Is  best 


considered  In  a cylindrical  coordinate  frame,  which  Is  fixed  In  space  and  It, 
therefore,  represents  the  Inertial  frame  of  reference.  The  classical  differential 
equations  of  motion  are  written  as 


■ F /m 

X 


r - r9  - F /m 
r 


r9  + 2r9  - F^/m 


where  (F^,F^,Fg)  are  components  of  the  net  applied  forces  on  the  roller  and  m 
Is  the  mass  of  the  roller. 


For  a given  position  of  the  roller  mass  center  and  all  translational  and  angular 
velocities,  the  applied  forces  are  calculated  from  the  various  Interactions  In 
the  bearing,  such  as: 


• Roller-race  Interaction 


e Roller  end  — race  flange  Interaction 

• Roller-cage  Interaction 

• Lubricant  drag,  churning  or  other  Interactions. 


The  mechanics  of  the  roller-cage  Interaction  have  been  formulated  In  (reference 
1) . The  details  of  roller-race  contact  and  the  roller  end  - race  flange  Interaction 


r 
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will  be  considered  In  Che  next  two  sections  of  this  report  and  other  possible 
Interactions  will  be  undertaken  as  a part  of  future  development. 

Once  Che  applied  forces  are  known,  Che  differential  equations  (II. 1)  of  motion 
can  be  readily  Integrated  to  obtain  Che  time  dependent  motion  of  Che  roller  mass 
center. 


Rotational  Motion  of  the  Roller 

If  the  roller  Is  considered  to  rotate  about  a fixed  point,  0,  with  angular 
velocity  (i)  and  the  triad,  formed  by  unit  vectors  3^,  ic  along  Che  principal 
axes  of  Inertia  at  0,  rotates  with  an  angular  velocity  then  the  classical 
differential  equation  of  motion  Is  written  as  (reference  7). 


■ 


X h 


- d 


(II. 2) 


where  Is  Che  roller  angular  momentum,  ^ 

ill  . T ^ 1 t ^2  in 

TF  “ ^1  dT  ^ ^ ^ 


Is  the  applied  moment  vector  and 

iT 


(11.3) 


where  Ij,,  1 2$  and  I^  are  the  principal  moments  of  Inertia  of  the  roller  . 

Some  straightforward  algebraic  manipulation  can  show  that  equations  (II. 2)  and 
(II. 3)  can  be  combined  to  give  the  following  component  equations: 


1,(1),  ~ 

" G, 

11  223  332 

1 

l2<i»2  - 

- S 

(II. A) 

I-<It«  “ I,u),J5»  + 

" G, 

33  112  221 

3 

In  case  of  a "perfectly”  cylindrical  roller,  the  transverse  moment  of  Inertia 

-*  I^  will  be  equal.  However,  this  symmetry  does  not  simplify  the  situation 
to  any  great  extent  since  the  polar  axis  has  to  be  fixed  In  the  roller.  It  will, 
therefore,  be  convenient  to  consider  the  roller  angular  motion  In  a roller  fixed 
frame  and  preserve  all  generalities.  A roller  fixed  frame  will  essentially  mean 
that  Che  triad  of  unit  vectors  1,  J,  ic,  of  equation  (II. 3)  Is  fixed  In  the  roller. 
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Thus,  n * u)  and  equations  (II. 4)  reduce  to  the  classical  Euler  equations  of 


motion 


“ °1 

12^2  - (I3  - - G2 


hS  • ^^1  - ^2>‘"l“2 


(II. 5) 


Clearly,  the  components  of  the  moment  vector  G will  now  be  In  the  body  fixed 
reference  frame. 

The  applied  moment  vector  G Is  once  again  determined  from  the  various  Inter- 
actions mentioned  In  the  case  of  roller  mass  center  motion.  It  Is  easy  to 
understand  that  while  the  normal  contact  forces  at  each  Interaction  dominate  the 
mass  center  motion,  the  angular  motion  Is  dependent  on  the  moments  which  result 
from  the  tractive  forces  at  the  various  interactions. 

Since  the  applied  moments  are  conveniently  expressed  in  the  Inertial  coordinate 
frame.  It  will  be  necessary  to  define  a law  of  transformation  between  the  inertial 
and  body  fixed  frames.  Although  such  transformations  have  been  defined  In 
(reference  1)  in  case  of  the  cage  motion  and  In  (reference  2)  for  ball  motion, 
they  will  be  repeated  here  for  completeness. 

Coordinate  Transformations 

As  shown  in  Figure  (II-2)  let  (X,Y,Z)  be  the  inertial  frame  and  (xb.yb.Zb^  denote 
the  body  fixed  frame.  The  required  transformation  is  obtained  in  terms  of  three 
angles  n,  3,  and  y by  the  following  three  successive  rotations: 


(1)  nt;  (11)  sr  (ill) 


where  I,J,K,l,j,k,  etc.  are  unit  vectors  along  the  corresponding  X,Y,Z,  and  x,y,z 


The  relevant  transformation  from  (X,Y,Z)  to  (^b’^b’S^  frame  may  be  written  as 


S ■ [Tjjj(n,B,Y)  1 R 


(II. 6) 
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where 


[T^^,(n,8,Y)] 


cosB  COSY 
-cosB  sinY 
sinB 


cosn  sinY 
+8 Inn  slnB  cosy 

cosn  COSY 
-slnn  sinB  sinY 

-Sinn  cosB 


Sinn  sinY. 
-cosn  sinB  cosy 

Sinn  cosy 
+co8n  sinB  sinY 

cosn  cosB 


(II. 7) 


It  may  be  noted  that  the  transformation  angles  n,  6,  and  Yt  are  quite  similar 
to  the  conventional  Euler  angles.  It  can  also  be  readily  shown  that  the  above 
transformation  matrix  is  orthogonal  and  hence  the  inverse  is  readily  determined 
by  the  transpose. 


In  order  to  complete  the  formulation  it  is  now  necessary  to  establish  the  rela- 
tionship between  the  angular  velocity  and  the  time  rate  of  the  transformation 
angles  n,  B,  and  y.  Let  w be  the  angular  vector  in  the  body  fixed  reference 
frame  (^>yjj»^^  ^^b’^b’^b^'  accordance  with  the  law  of  transformation 

defined  earlier,  the  incremental  vector  rotation  6n  is  written  as 


6n  - Anl  + ASJ"  + 


(II. 8) 


where 


T*  ■ SinY  ijj  + COSY  Jjj 

I ■ cos 3 t"  sinS 

■ cosB  (cosy  - siny  + sins 


(II. 9) 


A combination  of  the  above  two  equations  results  in  the  expression  for  re- 
quired angular  velocities 

6n  ■ (AncosScosY  + A68inY)t[j  + (-AncosgslnY  + AecosY)Jjj  + 

(AnslnB  +Ay)^ 
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b 

“1 

cosB  COSY  slnY  0 

-cosB  slnY  COSY  0 

b 

“3 

slnB  0 1 

^ J 

(II. 10) 


It  may  be  noted  that  the  matrix  In  Equation  (II. 10)  Is  not  orthogonal. 
However,  the  computation  of  the  Inverse  Is  rather  straightforward  and  some 
algebraic  manipulation  will  show  that 
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cosy 
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^ (II. 11) 
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Y 

V.  ^ 

-tanB  COSY 

tanB  slnY 

1 

1 

tt)  b 

3 

1 ^ 

The  accelerations  may  be  written  by  a straightforward  differentiation  of  the 
above  equation 
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where 
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and  Ic  may  be  shovm  that 


mi 

3t 


-Y  slny  sec8 
+5  COSY  sec3  tanS 


-Y  COSY  secS 
-8  sinY  secS  tanS 


Y cosy 


Y slny 


0 


siny  tanS  ^ 
COSY  sec^8 


Y 


COSY  tanB  2 
+3  siny  sec  S 


(11.14) 


For  prescribed  initial  conditions  the  above  formulation  may  be  readily  used 
to  obtain  the  real  time  simulation  of  the  ball  motion.  The  initial  conditions 
will  consist  of  ball  position  and  velocities.  Thus,  the  mass  center  position, 

(x,r,8),  mass  center  velocity,  (x,r,^),  angular  position,  (n,3,Y)  are  prescribed 

• • • 

at  any  instant  of  time.  Sometimes  w may  be  known  instead  of  (n.S.y),  in 
such  a case  equation  (II. 11)  will  provide  the  necessary  relationship.  With 
these  given  position  and  velocity  components  the  force  and  moment  vectors, 

F and  G respectively,  are  computed  from  the  models  for  the  various  inter- 
actions. Both  these  vectors  will  generally  be  determined  in  the  inertial 
frame.  The  mass  center  accelerations  are  given  by  equation  (II. 1).  For  the 
angular  acceleration  the  moment  vector  will  have  to  be  transformed  in 
the  roller  fixed  frame  using  equation  (II. 7)  and  (II. 8).  This  transformed 
moment  vector  can  then  be  substituted,  along  with  the  known  angular 
velocities,  in  equations  (II. 4)  and  oj  can  be  determined.  Equation  (11.12) 
is  now  used  to  compute  n,8  and  y.  All  the  acceleration  components  are 
finally  integrated  numerically  to  obtain  the  required  time  dependent 
motion. 
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SECTION  III 


ROLLER-RACE  INTERACTIONS 

The  normal  and  tractive  forces  at  the  roller-race  contact  constitute  the  domi- 
nant forces  on  the  roller.  The  Interaction  is  very  similar  to  that  in  case  of 
a ball-race  contact  (reference  2) . It  is  first  necessary  to  compute  the  geo- 
metrical interaction  between  the  roller  and  a raceway.  This  geometrical  inter- 
action will  define  the  elastic  deflection,  if  a contact  exists,  which  in  turn 
will  give  the  normal  contact  load,  using  a given  load-deflection  relationship. 
Once  the  contact  geometry  and  load  are  determined,  the  local  slip  velocities  in 
the  contact  zone  may  be  determined  from  the  prescribed  velocities  of  the  roller 
and  race.  A suitable  traction  model  for  a given  lubricant  and  operating  condi- 
tions is  then  used  to  compute  the  tractive  force.  A dlsciisslon  of  all  these 
factors  is  the  objective  of  this  section. 


Geometrical  Considerations 

Figure  (III-l)  schematically  describes  an  exaggerated  view  of  the  geometrical 

interactions  between  the  roller  anC  a race.  For  simplicity,  the  diagram  is 

drawn  in  the  x-z  plane  and  y axis  is  normal  to  the  plane  of  the  diagram.  The 

position  vector  r^  locates  the  mass  center  of  the  race,  R^,  relative  to  the 

inertial  frame  (X,Y,Z).  Vector  r locates  the  race  geometric  center,  R^,  rela- 

tive  to  the  mass  center  and  vector  r locates  the  center  of  race  land  R„ , 

gc 

relative  to  the  race  geometric  center,  R^.  Similarly  the  mass  center  of  the  roller, 

is  located  by  r^  relative  to  the  inertial  frame  and  r,  locates  the  roller  geo- 

metric  center  relative  to  the  ma-s  center.  The  vectors  r and  r.  will  generally 

rg  bg  * ^ 

be  prescribed  respectively  in  the  race  fixed  frame  (x,y,z)  and  the  roller  fixed  frame 
(x,y,z}.  Also,  the  relevant  transformations  will  be: 

Transformation  from  inertial  to  race  frame:  [T.  (n  ,6  ,y  )] 

ir  r’  r* 'r 

Transformation  from  inertial  to  roller  frame:  VI 

In  addition  to  the  above  it  will  be  convenient  to  define  an  azimuth  frame 
(x  ,y  ,z  ) such  that  z is  parallel  to  the  radial  component  of  r.  , x is  parallel 
to  the  inertial  X axis  and  y^  is  determined  by  the  right-hand  screw  rule.  Since 
the  diagram  in  Figure  (III-l)  is  drawn  in  the  x-z  plane,  this  coordinate  frame  is 

Ill-l 


Figure  III-l.  Exaggerated  View  of  Roller-Race  Interaction 
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not  shown  In  Che  figure.  However,  the  transformation  will  depend  only  on  the 
race  azimuth  angle  Y,  and  it  may  be  defined  as 


1 

i 


I 
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Transformation  from  inertial  to  roller  aximuth  frame:  [T^^(Y,0,0)} 

In  order  to  make  the  entire  analysis  applicable  to  both  the  outer  and  inner  race 
contacts,  it  will  be  convenient  to  introduce  a transformation  from  the  azimuth 
to  contact  frame. 

Transformation  from  azimuth  to  contact  frame:  (T^^(0,a,0)] 
where  a > 0 fcr  outer  race  and  for  inner  race. 

It  may  be  noted  chat  for  most  practical  problems  the  mass  and  geometric  centers 
may  be  coincident  and  this  could  simplify  the  geometry  Co  some  extent.  The 
purpose  of  defining  Che  geometric  center  different  from  the  mass  center  is  just 
to  preserve  Che  generality  and  capability  of  simulating  some  advanced  config- 
urations in  the  future. 


Now  the  interaction  between  the  roller  and  Che  race  will  be  determined  by 
locating  the  geometric  center  of  the  roller  Bg,  relative  to  the  race  lend  center. 
Let  us  denote  this  vector  by  r^^^  and  it  will  be  given  by  the  equation 


+ [T 
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♦ <'rg 
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+ r 

gc 


(III.l) 


It  may  be  noted  that  everything  on  the  right-hand  side  of  Che  above  equation 
is  known  except  for  Che  vector  r . It  will  be  necessary  to  determine  a race 
azimuth  angle  measured  in  the  race  fixed  coordinate  frame  such  that  B^,  R^, 
and  Rg,  all  lie  in  one  plane.  The  vector  locating  B^,  relative  to  R^, 

written  as 


•♦•1  -*■ 

■ r.  + [T,.  Jr. 
brg  b ib  bg 


(III. 2) 


When  ^^j.g  Is  1°  inertial  frame,  an  angle  may  be  defined  such  Chat  it 

denotes  Che  angle  between  Che  inertial  z axis  and  Che  plane  containing  vector 

r.  and  a line  passing  through  R-  and  parallel  Co  Che  inertial  X axis, 
org  v» 
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Symbolically 


r-r, 


y ■ arctan 


brg. 


brg- 


(III. 3a) 


If  Che  race  is  assumed  to  be  symmetric  about  the  x axis,  then  the  race  fixed 

coordinate  frame  may  not  be  necessary  and  $'  - f.  However,  if  the  generalities 

mentioned  above  are  preserved,  then  , is  determined  in  the  race  fixed  frame 

such  that  r and  r,  lie  in  the  same  plane.  Now  if  \ is  the  radius  of  the 
gc  brg 

circle  representing  the  locus  of  the  race  land  centers,  then 
{0,  -<J.sin4>' , Rcos  *} 


-►r 

r 

gc 


In  order  Chat  r^^  and  are  coplaner,  the  y-component  of  the  following  vector 

must  vanish 
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or  » arctan 


23 

\'^22, 


where  T22  and  relevant  components  of  the  matrix  obtaired  by  the 

product  [T('i',0,0)] 

Once  $ is  known,  may  be  determined  in  the  inertial  frame  using  equation 
(IH.l)  and  it  may  be  referenced  in  the  race  azimuth  frame  as 


br 


[T(«’,0,0)]r, 
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and  in  Che  contact  frame 
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where  the  matrix  [T  ] has  already  been  defined  above, 
ac 

The  computation  of  the  roller  azimuth  tj),  is  quite  straightforward  once  the 
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orbital  poalcion,  6,  locating  th«  rnllar  maaa  cantar  la  kno«m  from  praacrlbed 
r^  3 (x.r.O)  In  cylindrical  coordlnatas. 


li!  ■ 0 - arc  tan 


(III. 6) 


whara  r2  and  r ^ ara  tha  y and  r componants  of  tha  vactor  Ir^  . 


Onca  tha  gaomatrlc  cantar  of  tha  roller  la  located  relative  to  the  race  land 
center,  It  will  be  neceasary  to  examine  every  point  on  the  roller  for  poaalble 
Interaction.  Hanca,  the  Introduction  of  the  roller  geometry  becomea  easential 
at  this  point.  Figure  III-2  showa  this  geometry  for  a typical  roller  In  a 
cylindrical  roller  bearing.  The  Interaction  between  the  roller  and  race  will  be 
restricted  to  tha  central  flat  land  and  the  crowned  areas  at  this  point.  Later 
In  the  next  section  the  Interaction  between  the  roller  end  and  race  flange  will 
be  considered. 

Now  a point  of  the  roller  surface  Is  bast  described  In  Che  roller  frame  by  Che 
position  vector 
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m ,-R  sin.j 

R-dosiJ 


(III. 7) 


Where  i'  is  Che  angular  poalclon  about  the  X axis,  I Is  Che  length  coordinate 
along  the  X axis  and  R'ls  Che  effective  radius.  For  Che  geometry  shown  In 
Figure  III-2,  will  depend  on  4 

R-  . d.  . h < I < h 
R • 2*  2 - - 2 
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R'-  2 - -4r~ 

c 


^2  ^2 
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where  R^  la  the  crown  radius  as  shown  In  Figure  1II-2. 


The  vector  r can  be  expressed  In  the  contact  frame  as 

r - (T  1 (T  1 (T,  1 (T*  1 r^ 
ac  ‘ ra'  ‘ Ir'  lb 


(III. 8) 
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Figure  III-2.  Roller-Race  Contact  Geometry 
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For  the  point  to  be  In  the  contact  plane  the  y component  of  ^ must  vanish 
and  t-'  '3  condition  will  determine  Che  relevant  Since  the  rotation  of  the 
roll  about  Che  Y and  Z axes  will  be  in  general  small,  it  will  only  be  neces- 
sary Co  compute  ^ for  the  roller  center,  hence 

T 

♦ - arctan  (III.  9) 

22 

where  and  T22  corresponding  elements  in  a matrix  obtained  by  the 

product  of  matrices  shown  in  equation  (III. 8). 


The  deflection  at  any  point  on  the  roller  surface  is  now  given  as 


5 - (r^  + rf  ) (III. 10) 

br  3 


A negative  value  of  5 will  mean  no  contact. 


In  order  to  determine  the  effective  contact  length,  when  the  roller  is  crowned, 
Che  length  coordinate  I in  equation  (III, 10)  will  be  determined  by  standard 
bisection  techniques.  In  other  words,  Che  value  of  I for  5 ■ 0 will  determine 
the  contact  length.  Clearly,  these  will  be  two  such  roots  and  they  both  have 
to  be  determined  separately  since  in  general  the  roots  may  not  be  symmetric. 

For  the  purpose  of  computing  the  total  interaction,  it  will  be  convenient  to 
divide  the  contact  length  into  a finite  number  of  elements  as  shown  in  Figure 
III-2  and  compute  the  normal  and  tractive  loads  for  each  element.  A proper 
summation  will  then  give  the  total  required  load.  Thus,  the  analysis  presented 
below  for  Che  normal  and  tractive  load  will  apply  to  each  elemc  nC  in  the  contact 
zone. 


Normal  Contact  Load 


The  load  deflection  relation  for  a line  contact  is  not  uniquely  defined,  as  in 
case  of  point  contact.  The  deflection  can  be  computed  only  in  terms  of  strains 
corresponding  to  finite  lengths  in  Che  two  interacting  bodies.  (Reference  3} 
suggests  a relationship  of  the  form 
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where  the  contact  half  width 
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(III. 12) 


'^1’  ^1’  '^2’  ^2’  Poisson's  ratio  and  elastic  modulus  for  the  two  inter- 

acting bodies.  Ip  is  the  sum  of  the  principal  curvatures,  Q is  the  normal 
contact  load,  and  I is  the  contact  length. 

The  determination  of  lengths  d^  and  62  is  rather  arbitrary  and,  hence,  application 
of  (III. 11)  becomes  a little  difficult. 

In  view  of  this  difficulty,  Lundberg  (reference  4)  has  determined  an  empirical 
equation  based  on  actual  load  deflection  measurements  in  line  contact.  This 
relation  is 
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where  the  contact  half  width  b is  given  by  equation  (III. 12). 

It  is  clear  that  for  a given  Q,  the  computation  of  6 is  straightforward  from 
equation  (III. 13).  However,  when  6 is  prescribed,  some  type  of  iterative  pro- 
cedure will  be  necessary  to  order  to  compute  Q. 

Palmgren  (reference  5)  gives  a relation  which  is  free  of  the  above  difficulty. 


0.39 
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(111.14) 


In  terms  of  results,  the  above  two  relations  are  fairly  close.  Hence,  if 
equation  (III. 13)  is  to  be  used,  then  the  initial  guess  in  the  iterative  scheme 
can  be  derived  from  equation  (III. 14)  and  under  such  conditions  the  Iteration 
will  converge  rather  rapidly. 

Tractive  Forces  and  Moments 

For  the  purpose  of  computing  the  traction,  it  will  be  necessary  to  define  the 
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position  vector  locating  the  center  of  the  elementary  strip  In  the  contact  zone, 
relative  to  the  roller  mass  center  and  the  local  slip  velocities.  If  we  assume 
that  both  roller  and  race  deform  equally,  which  Is  quite  reasonable,  then  the 
required  position  vector  may  be  written  as 


- |-(R-6/2)sin(>  ( 

l(R-6/2)co3<t  J 


(III. 15) 


In  order  to  compute  the  local  slip  velocities,  let  the  race  have  a translation 
velocity  v^  and  an  angular  velocity  Likewise,  let  the  roller  have  a trans- 
lational velocity  Vj^  and  an  angular  velocity  Generally,  the  angular 

velocities  are  specified  in  the  body  fixed  frames  while  the  translational  veloc- 
ities are  measured  in  the  inertial  frame.  Furthermore,  the  roller  velocity  is 
conveniently  measured  in  the  cylindrical  coordinate  frame  in  terms  of  components 

X,  f,  and  9.  In  terms  of  these  components,  the  local  slip  velocity  u of  the 

s 

race  relative  to  the  roller  at  point  P (see  Figure  III-2)  Is  given  In  the  contact 
frame  in  terms  of  the  race  and  roller  velocities,  u^  and  u^^  respectively  as 
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where  the  vector  R , which  locates  the  center  of  elementary  strip  with  respect 
to  the  race  mass  center.  Is  computed  as 
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(III. 19) 


c c c 

The  components  u and  u in  the  contact  frame  will  result  in  traction,  u 

®1  ®2  c ®2 

being  along  the  rolling  direction  and  u^  being  normal  to  it. 

®1 


In  case  of  a roller  bearing  the  slip  will  be  primarily  in  the  y direction  and, 
hence,  the  tractive  force  will  act  along  the  y axis.  The  relevant  traction 
coefficient,  k,  will  be  determined  by  the  appropriate  traction  model,  which  will 
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be  described  In  Section  V,  for  the  present  It  mey  be  denoted  as; 


*2  ” 


where  the  subscript  2 is  used  to  denote  the  y direction. 


The  tractive  force,  F^,  Is  now  defined  as 


Fj.  - ic:} 


(III. 20) 


(III. 21) 


where  Q Is  the  already  known  normal  load  on  the  elementary  strip. 

Both  the  normal  and  tractive  forces,  acting  on  the  roller  can  be  combined  in 
a general  force  vector,  written  conveniently  in  the  contact  frame  as 


The  force  on  the  race  will  just  be  -F,  while  the  moments  on  the  roller  and 
race  are  defined  as 


S ■ "p  * ^ 


and  G «(R+r)xF 
r p rg 


All  vectors  must  be  transformed  to  appropriate  coordinate  frames.  It  is  con- 
venient to  compute  the  forces  in  the  initial  frame  and  the  moments  In  a body  fixed 
frame. 

It  should  be  remembered  that  all  the  above  forces  and  moments  will  be  on  a given 
elementary  strip  In  the  contact  zone  and  a summation  over  the  contact  zone  will 
give  the  total  force  and  moment  vectors. 


III-IO 


SECTION  IV 

ROLLER  END-RACE  FLANGE  INTERACTION 


In  a roller  bearing  subjected  to  a slight  mlsallgiment , the  loads  between  the 
roller  and  raceways  may  be  somewhat  nonsymmetric  and,  hence,  the  resulting 
tractions  develop  a moment  about  the  transverse  axis  of  the  roller.  Such  a 
moment  will  lead  to  angular  rotation  of  the  roller  about  its  transverse  axis 
and  this  type  of  motion,  which  Is  highly  undesirable.  Is  generally  called 
skewing.  It  is  quite  easy  to  understand  that  when  the  roller  skews,  the  corners 
of  the  roller  could  interaction  with  the  cage  and  the  flanges  on  the  race. 
Generally,  the  axial  clearance  between  the  roller  end  and  cage  is  quite  large 
compared  to  that  between  the  roller  end  and  race  flange.  Hence,  the  cage  inter- 
action can  be  neglected,  but  It  Is  essential  that  the  roller  end  and  race  flange 
Interaction  be  modelled  in  a realistic  fashion  in  order  to  treat  the  roller  skew- 
ing with  adequate  justice.  Unfortunately,  the  local  geometry  and  the  kinematics 
greatly  complicate  the  problem  and  the  development  of  a realistic  model  becomes 
extremely  difficult.  It  should  be  remembered  that  a realistic  model  should  not 
only  simulate  the  problem  correctly,  but  it  should  also  be  simple  enough  to  be 
rapidly  executable  of  a digital  computer,  since  computer  time  Is  a key  factor 
In  sophisticated  programs,  such  as  the  one  under  consideration. 

With  the  above  difficulties  in  view,  a first  order  interaction  model  for  the 
roller  end  and  race  flange  contact  is  derived  on  the  basis  of  the  following  two 
assumptions: 

(1)  If  the  roller  end  is  not  contacting  the  race  flange,  there  is  no 
contact  force.  In  other  words,  any  hydrodynamic  action  is  neglected. 

(2)  When  the  race  flange  is  contacting  the  roller  corner,  which  will  have 
a finite  radius  of  curvature,  the  contact  load  is  approximated  by  an 
effective  Hertzian  type  contact  model  and  the  tractive  force  is  deter- 
mined by  a prescribed  traction-slip  relation. 

It  is,  therefore,  clear  that  the  objective  is  to  treat  metal  to  metal  contact  only. 
There  will  basically  be  two  parts  to  such  a treatment,  e.g.,  the  contact  geometry 
and  the  load  deflection  relationship.  Both  of  these  parts  will  be  discussed  in 
this  section. 
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Geometrical  Consideration 


The  roller  geometric  center  was  located  with  respect  to  the  race  geometric  center 
In  the  preceedlng  center.  In  case  of, the  roller  corner  and  flange  interaction, 

t 

It  will  be  necessary  to  locate  the  center  of  curvature  of  the  roller  corner  with 
respect  to  a point  on  the  race  flange.  Thus,  an  additional  transformation  from 
the  race  azimuth  to  flange  frame  la  Introduced. 


Transformation  from  race  azimuth  to  race  flange  frame:  (T^^  (0,y,0))  where  the 
angle  y will  depend  on  the  particular  flange  under  consideration  as  shown.  In  an 
exaggerated  fashion,  in  Figure  lV-1. 
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It  should  also  be  noted  chat  the  race  azimuth  angle  will  now  be  defined  by  the 


relevant  center  of  curvature  of  the  roller  corner  and  not  by  the  center  of  roller. 


Let  us  consider  a point,  0,  on  the  locus  of  the  center  of  curvatures  of  the  roller 
corner,  as  shown  In  Figure  IV-2.  It  Is  clear  that  the  locus  will  be  a circle  and 
the  position  vector  locating  0,  with  respect  to  the  roller  geometric  center  in 
roller  frame,  is  given  as  1 
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where  is  just  angular  coordinate  about  the  X axis  in  the  roller  fixed  frame. 


It  will  be  necessary  to  write  r^  in  the  race  frame,  and  ultimately  in  the  flange 


frame,  and,  hence,  we  recall  the  transformation  from  the  roller  to  race  frame  from 
the  preceeding  section 
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Also,  a vector  locating  the  roller  geometric  center  with  respect  to  the 


race  geometric  center  can  be  written  In  terms  of  the  various  vectors  Introduced 
In  Figure  III-l  in  the  preceeding  section. 
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Figure  IV-2.  Exaggerated  Viev  of  Roller  Comer  and  Race  Flange  Interact  ion 
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Nov  for  any  selected  value  of  1 In  (IV. 1),  the  race  azimuth  angle  ip  Is  defined, 
such  that  the  vector 

r ■ r + r. 
er  e brg 

when  written  In  the  race  azimuth  frame  has  only  X and  Z components.  Hence, 

r ■ r,  + {T,  } r^*^ 
er  brg  br  e 


and  r 


(IV. 4) 


where 


a r . a 

r » r and  r 
er^^  er  ^ er^ 


. r .2  . , .2 

(r  ) + (r  ) 

er2  er^ 


It  should  be  remembered  the  r locates  the  roller  corner  curvature  center  with 

er 

respect  to  the  race  geometric  center.  If  we  define  the  point  O'  as  a reference 
point  on  the  flange,  then  the  roller  curvature  center  can  be  referenced  relative 
to  O'  as 


— ►a  — ►a 

r , * r 

ef  er 


■s 


(IV. 5) 


where  S Is  the  race  half  width  and  d Is  the  radius  of  the  race  land. 


When  the  above  Is  written  In  the  flange  frame 


(IV. 6) 


Clearly,  the  point  of  maximum  Interaction  will  be  determined  by  maximizing  the 
Z component  of  fact,  the  Z component  will  define  the  elastic  deflection 

If  It  Is  positive,  and  greater  than  the  comer  radius, 

It  may  be  recalled  that  all  of  the  above  Is  carried  out  for  a preselected  value 


1 


r 


of  ^ in  equation  (IV. 1).  Thus,  (r  -)  is  maximized  by  varying  4i . Hence, 

et  3 


After  carrying  out  some  routine  algebra,  it  may  be  shown  that  the  above  equation 
may  be  written  as 


^rf„  f-  '^br_  ^ 


rf  r 3 

+ — I I (r^  + ’^br  ^ " ^br  ^ ^ 

r^r  ii=2  ’"’"i  - ‘’’'il  “i2  ^’'13 


X R cos,^  - R sin^)  > - 0 


(IV. 7) 


The  above  is  solved  for  ? by  standard  bisection.  Note  that  there  will  be  two 
possible  roots,  one  denoting  the  maximum  and  the  other  defining  a minimum.  The 
interest  here  is,  of  course,  in  the  maximum. 


As  discussed  above,  the  elastic  deflection  will  be  defined  by  the  (r^^  ) when 
it  is  positive.  ^ 


5 = (r  c ) ~ H when  >0 

ef  - ce 

3 max 


(IV. 8) 


If  6 is  negative,  then  there  will  be  no  interaction  between  the  roller  corner 
and  the  race  flange. 


Normal  Contact  Load 

If  the  contact  between  the  roller  corner  and  race  flange  exists,  then  the  contact 
load  can  be  determined  by  using  the  classical  point  contact  Hertzian  solution 
described  earlier  (reference  2),  Clearly,  the  deflection  5 is  computed  as 
described  above,  the  material  properties  are  prescribed  and  the  only  additional 
parameters  needed  are  the  effective  radii  of  curvature  of  the  interacting  bodies. 
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which  may  be  defined  as  follows. 


Tr,  ■ « 

Race  flange  ^ ^ 


Roller  corner 


1*^1  - \e 

[^2  “ ^^ce  ■*■  R)/‘^oseg 


where  9^  is  the  skew  angle.  This  cingle  is  Included  in  the  transformation 
used  above  and  in  the  preceeding  section. 

It  should  be  noted  that  when  6 is  zero,  the  above  model  will  not  be  valid, 

s 

since  the  entire  end  face  of  the  roller  will  interact  (for  a flat  end)  and  the 
contact  is  no  longer  Hertzian.  Such  a contact  will  only  exists  when  there  is  a 
dominant  thrust  load  on  the  bearing,  which  is  never  the  case  in  cylindrical  roller 
bearing. 


An  alternate  way  to  compute  the  normal  load  will  be  similar  to  the  one  described 
for  roller-cage  and  ball-cage  interactions  in  (reference  1).  A certain  stress 
level  for  the  computed  deflection  is  prescribed  and  the  contact  area  is  just  com- 
puted by  the  geometry  of  interaction. 


Tractive  Forces  and  Moments 

The  computation  of  tractive  force  will  be  very  similar  to  the  procedure  described 
for  the  roller-race  interaction  presented  in  the  preceeding  section.  The  contact 
point  is  located  relative  to  the  roller  and  race  mass  centers  and  the  slip 
velocities  are  computed.  So,  if  r^  defines  the  contact  point  relative  to  the 
roller  mass  center,  and  R^  defines  it  relative  to  the  race  mass  center,  it  will 
be  straightforward  to  show  that 


ce\ 


““  * C'  * 'rt 


(IV. 9) 
(IV. 10) 


where  [T  1 is  the  transformation  from  the  race  frame  to  the  race  azimuth 
ra' 

frame. 
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The  local  velocities  at  the  contact  point  are  clearly 


% ” '"b  “b  ""p 


u • V + u)  X R 


r r r p 

Effective  slip  u - (u  - u.  ) 
s r2  *^2 

and  the  traction  coefficient  k = 


'^(U3) 


U 

s 


(IV. 11) 
(IV. 12) 


The  total  force  acting  on  the  roller  is  now  defined  in  the  flange  frame  as 


where  Q is  the  normal  contact  load,  discussed  above. 


(IV. 13) 


1 of  the 

race  will  just  be 

“ - K 

(IV. 14) 

moments 

on  the  roller  and  raceway  will  be 

- -+• 

'bg’  X ^ 

(IV. 15) 

^ 

=■  (R  + 
P 

Xfg  X ffc 

(IV. 16) 

It  must  be  remembered  that  all  vectors  must  be  transformed  in  appropriate 
coordinate  frame.  Generally,  the  forces  will  be  computed  in  the  inertial  frame 
and  the  moments  in  the  body  fixed  frame  of  reference. 
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SECTION  V 


ELASTOHYDRODYNAMIC  TRACTION  MODELS 

Several  approaches  to  a realistic  model  of  the  traction  slip  relationship  of 
a given  lubricant  may  be  considered.  The  simplest  approach  will  be  to  test 
the  lubricant  at  the  prescribed  operating  condition  and  derive  a hypothetical 
model  by  curve  fitting  the  experimental  data.  A somewhat  more  realistic 
approach  will  be  to  perform  a two-fold  elastohydrodynamic  analysis  where  the 
first  part  will  determine  the  nominal  film  thickness  and  the  second  part  will 
evaluate  the  effective  values  of  the  coefficients  in  the  pressure-temperature- 
viscosity  relation  by  curve  fitting  the  experimental  traction  data.  A discussion 
of  these  different  approaches  is  the  primary  objective  of  this  section. 

A Hypothetical  Model 

Most  of  the  available  experimental  data  shows  a general  trend  demonstrated  in 
Figure  II-l.  Traction  initially  increases  with  Increasing  slip,  reaches  a 
maximum  value  at  some  slip  velocity,  and  finally  reduces  to  some  asymptotic 
value  as  the  slip  continues  to  Increase  further.  Such  a qualitative  trend 
may  be  simulated  by  an  expression  of  the  form  (reference  8) 


< - + 1^2  u)  exp  u)  + 


(V.l) 


where  < is  the  traction  coefficient  and  u is  the  slip  velocity. 


In  general,  the  coefficients  Pjt  and  will  be  a function  of  the 

operating  conditions  and  the  lubricant  properties.  However,  if  it  is 
postulated  that  the  model  is  valid  for  a given  lubricant  at  prescribed  oper- 
ating conditions,  then  the  values  of  the  coefficients  may  be  derived  from 
the  various  parameters  shown  in  Figure  V-l,  which  can  be  expressed  as 


u • 0, 


< ■ K. 


U ^ 

m 


3<  . 

< 1 ■ 0 

m 3u 


(V.2a) 

(V.2b) 


(V.2c) 


Traction  Coefficient 


r 

i 

[ Combining  (V.l),  (V.2a),  and  (V.2b)  gives 

(V.3) 
(V.4) 

Also,  from  equations  (V.l)  and  (V.2c) 


i 

I 


-'P3  + 'l>2  %)  + '<'2  ■ 0 (V.5) 

and 

(i',  + (p-  u ) exp  (-Ip  u ) + Ip,  - K i 

X ^ to  j n 4 IS  i 

which  form  the  two  simultaneous  equations  for  determining  1P2  '^2’ 

Eliminating  in  equations  (V.5)  will  give  a nonlinear  equation  in 

S 

i 

(V.6)  I 

which  is  solved  by  standard  bisection  methods.  Hence,  all  the  four  coefficients 
can  be  computed  from  the  prescribed  conditions  shown  in  Figure  V-1. 

It  is  clear  that  the  above  model  is  very  simple  in  that  it  does  not  allow 
any  variation  in  the  parameter  as  a function  of  local  variations  in  rolling 
speed,  Hertz  stresses,  etc.  In  order  to  allow  for  these  variations,  a some- 
what more  sophisticated  elastohydrodynamic  model  will  be  required.  The  first 
step  in  such  a model  is  the  computation  of  the  nominal  film  thickness,  which 
will  be  discussed  below. 

Nominal  Film  Thickness  Computation 

The  computation  of  film  thickness  is  one  of  the  most  fundamental  steps  in  j 

simulating  the  behavior  of  a concentrated  elastohydrodynamic  contact  in 
rolling  bearings.  Most  of  the  available  analysis  first  determines  the  film 
thickness  under  isothermal  conditions  and  the  computed  value  is  later  modified 
by  a thermal  reduction  factor  in  order  to  accommodate  the  thermal  effects. 

Some  of  the  formulae  and  computational  procedures  are  discussed  below. 
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Grubln  Formula;  This  is  perhaps  one  of  the  first  relations  which  appeared 
in  the  literature  (reference  9 ) for  the  computation  of  the  isothermal  film 
thickness  h,  in  a line  contact  configuration 


h 

R 


1.95 


(GU) 

W 


8/11 

iTir 


where 


R 


^ --p  . the  effective  radius  in  rolling  direction 

R1-R2 


(V.7) 


G » aE',  the  elasticity  parameter 


U 


the  speed  parameter 


U 


i '^2^’ 


the  rolling  speed 


W = . the  load  parameter 

V,  E,  a,  and  w are  respectively  the  Poisson's  ratio,  elastic  modulus, 
inlet  viscosity,  pressure  viscosity  coefficient,  and  applied  load  per  unit 
length. 


Dowson-Higginson's  Formula;  A similar  relation  is  put  forth  by  Dowson  and 
Higginson  (reference  9) 


h 

R 


0.60y0.70 


(V.8) 
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The  various  parameters  are  the  same  as  these  described  in  the  case  of  the 
Grubln's  formula. 

Thermal  Reduction  Factor;  The  isothermal  film  thickness  computed  by  any  of 
the  above  expressions  must  be  modified  by  a thermal  reduction  factor  in 
order  to  allow  for  appropriate  thermal  effects.  Cheng  (reference  10)  has 
summarized  the  results  of  a thermal  analysis  and  the  required  thermal 
reduction  factor  can  be  computed  by  interpolating  the  results.  Primarily 
for  completeness,  these  results  are  reproduced  here  from  reference  10. 


The  thermal  reduction  factor,  is  expressed  In  terms  of  the  following 
parameters 

“o  2M2i_ 


2 K-  T 
f o 


K,  T 
f o 


a*  “ Y . 10^  . a 


V2-V1 


and  — ^ , the  dimensionless  Hertzian  pressure. 

£. 


Kg,  6,  and  T^  are  the  thermal  conductivity,  temperature  viscosity  coefficient, 
and  inlet  temperature,  respectively.  All  other  parameters  are  the  same  as 
discussed  earlier. 


The  solutions  of  ij)_  indicate  that  the  most  influential  parameters  are  the 

^ ^Hz 

viscous  heating  parameter  Q^.  The  other  four  parameters,  a*,  6 , -gr.  and  s 

all  play  secondary  roles  in  governing  the  thermal  reduction  factor.  In 
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particular,  the  influence  of  the  slip  ratio  is  very  mild.  A review  of  these 
curves  shows  that  within  the  range  of  the  parameters  investigated,  the  thermal 
reduction  factor  can  be  represented  by 


I 


(1  - 0.1  s)(l  - 


(V.9) 


where  f,  and  are  numerical  functions  of  Q , a*,  8',  and  are  tabulated  in 
L z m 


Tables  V-1  and  V-2, 


Side  Leakage  Factor;  Cheng  (reference  10)  has  also  presentta  a side  leakage 
factor,  <ti^,  to  modify  the  line  contact  film  thickness  in  the  case  of  ellipti- 
cal contacts  which  may  be  relevant  of  crowned  rollers.  The  dimensionless 
film  thickness  for  an  elliptic  contact  is  expressed  as 


h r p aU  T n,  rPu,“|  n. 


(V.IO) 


where  C,  n^^,  and  n^  are  numerical  functions  of  a/b.  These  values  are  shown 


in  Table  V-3. 


Assuming  that  the  film  thickness  at  k - a/b  = 5 approaches  that  of  a line 
contact,  then  the  side  leakage  reduction  factor  becomes 


C, 


2 " ^”2^k=5] 


(V.ll) 


The  effective  film  thickness  can  now  be  expressed  as 


H * Vi  6— 

iso  ^T  S 
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TABLE  V-1 
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VARIATIONS  OF  WITH  RESPECT  TO 


* S'  . » 
a , — , and  Q 

a " 


■ a • -y  • 10^,  B'/a*  - — . i 
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15.71 

1 

15.71 

15.71 

23.564 

0.35 

0.5 

0.75 

0.35 

1.0 

1.0 

1.0 

1.0 
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.98 

.99 

.99 
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.93 

.9 
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.9 
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.86 

.81 

.86 

.8 

.78 

.66 

.77 

.^2 

.65 

.56 

.67 

.62 

.53 

.44 

.57 

.47 

.37 

.3 

.44 

.36 

.27 

.22 

.34 



.97  .97 

.95  .94 

.92  .90 

.87  .84 

.81  .76 

.7  .62 

.6  .50 
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31.419  31.419  31.419 
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.98 

94 
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.95 
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00 

84 

.90 

00 

.82 

76 

.84 
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.74 

62 

.75 

.66 

.59 

50 

.66 

.56 

.46 

36 

.56 

.45 

.35 

24 

.43 

.31 

.23 

17 

.34 

.21 
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TABLE  V-2 


VARIATION  OF  WITH  RESPECT  TO 
a*.  S' /a*,  and 
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15.71 

15.71 

15.71 

23.56 

23.56 

23.56 

31.42 

31.42 

0.35 

0.5 

0.75 

0.35 

0.5 

0.75 

0.35 

0.5 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

3.6 

3.95 

0.0 

3,78 

4.34 

0.0 

8.3 

0.0 

12.8 

12.4 

0.0 

11.4 

8.2 

4.1 

12.4 

47.5 

48.8 

40.7 

26.5 

17.6 

30.1 

12.9 

13.75 

107.8 

75.6 

51.4 

52.9 

26.5 

25.6 

25.6 

20.9 

TABLE  V-3 

NUMERICAL  FUNCTIONS  OF  C,  AND  n2 


k “a/b 

C 

"l 

"2 

5 

1.625 

0.74 

-0.22 

2 

1.56 

0.736 

-0.209 

1 

1.415 

0.725 

-0.174 

0.5 

1.132 

0.688 

-0.066 

The  Turbo-33  Oil  Model 


1 


Based  on  the  experimental  data  (reference  11)  on  the  Shell  turbo- 33  type  oil, 
a traction  model  has  been  presented  in  reference  10.  It  is  shown  that  the 
traction  coefficient  <2>  rolling  direction  is  dependent  on  three  param- 

eters 


8K^ 


ap 


H 


(V.13) 


where  is  the  lubricant  viscosity  at  an  inlet  temperature  of  T^;  is 

the  thermal  conductivity  of  the  lubricant;  h is  the  film  thickness;  U2  is  the 
slip  velocity  in  the  rolling  direction;  Pjj  is  the  maximum  Hertz  pressure;  and 
the  coefficients  a and  8^  are  obtained  from  the  viscosity-pressure-temperature 
model  of  the  form 


u(p,T)  • exp  [ap  - 8j^(T  - T^)[l 


(V.14) 


The  relationship  between  G2,  G^,  and  the  traction  coefficient  <2  is 
represented  graphically  (reference  10)  by  a series  of  graphs  and  this  data 
is  basically  stored  in  a computer  data  file  for  real  applications.  The 
film  thickness  h is  computed  by  the  methods  described  above. 


With  the  assumptions  of  a narrow  contact  ellipse  this  model  is  directly 

applicable  to  a crowned  roller  contact.  In  equation  (V.13),  P is  replaced  by 

H 

the  pressure  as  a function  of  S. 


and  U2  is  clearly  a function  of  C,  the  coordinate  along  the  major  axis  of 
the  contact  ellipse  (see  Figure  V-2).  Thus,K^  is  obtained  as  a function 
of  C. 
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The  5P4E  Polyphenyl  Ether  Model 

Smith,  et  al  (reference  12)  have  proposed  a traction  model  based  on  the  experi- 
mental data  obtained  by  them  for  polyphenyl  ether.  It  is  shown  that  the  shear 
stress  in  the  lubricant  is  best  correlated  by  an  expression 


sin  ~ [i//  exp  (a*  p/2)]  exp  (a*  p/2) 

h tj/  /l  + exp(a*P) 


(V.15) 


where 


p is  the  local  pressure  obtained  by  Hertzian  pressure  distribution,  (5,?)  is 
a coordinate  system  in  the  contact  ellipse  along  the  major  and  minor  axes  (see 
Figure  II-2) , and  other  parameters  are  the  same  as  defined  earlier.  The 


three  empirical  constants,  u 


B , denote  the  coefficients  in  a modified 


pressure-temperature-viscosity  relation 


exp  {a*p  + B*  (T  - T^)} 


(V.16) 


and  the  value  of  these  constants  are  determined  by  correlating  the  experi- 
mental data. 

With  the  assumption  of  narrow  ellipse  equation  (V.15)  may  be  integrated  with 
respect  to  C to  obtain  a tractive  force  per  unit  length  along  the  major  axis 
C.  If  this  force  is  divided  by  the  normal  force  per  unit  length,  a traction 
coefficient  rolling  direction  determined.  Some  straightforward 

algebraic  manipulation  will  show  that 
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iJiCC.C)  ■ 4i  exp 


and 


be  - b 


/l  - 


Generally  the  constants  a and  B are  independent  of  rolling  speed,  while 
u^*  is  allowed  to  vary  slightly  with  the  rolling  speed  in  order  to  allow  for 
the  variation  in  traction  as  a result  of  short-time  effects  on  lubricant 
viscosity  or  the  inlet  zone  heating  effects  (reference  12).  it  is  found 

(reference  12)  that  the  values  for  the  three  coefficients  which  fit  the  experi- 
mental data  are 


* —3  2 

U ■ 1.4  X 10  Ibf-sec/in  , for  rolling  speed,  U < 900  in/sec 


1.4  X 10“^  exp  1-0.30111 


f-L- 

l900 


Ibf-sec  ~ ^ ___  , , 

T — , U > 900  in/sec 

in 


a*  * 3.77  X 10"^  in^/lbf 


and 


0.046  1/OR 
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It  should  be  noted  that  relation  (V.16)  is  valid  only  for  the  computation 
of  traction.  For  determining  the  film  thickness  a relationship  of  the  form 


(V.18) 


is  used,  in  which  a,  0,  and  y are  the  conventional  pressure-viscosity  temperature- 
viscosity,  and  pressure-temperature-viscosity  coeffielcnts.  For  5P4E  poly- 
phenyl ether,  these  constants  have  the  following  values. 

T - 630  °R  U - 4.748P  Ibf-sec/in^ 

o c 

a - 1.1597  X 10““^  in^/lbf  6 - 6.5042  °R 

Y - 0.18305  °R  Ibf/ln^ 


Using  the  above  values,  the  inlet  viscosity  at  any  temperature  can  be  deter- 
mined. An  additional  property  required  in  computing  the  thermal  reduction 
factor  is  the  thermal  conductivity  which  is  found  to  be  1.205  x 10  lbf/sec/°F. 


The  MIL-L-7S08  Model 


Similar  to  the  polyphenyl  ether  model,  Walowit  and  Smith  (reference  13)  have 
investigated  the  behavior  of  the  MIL-L-7808  oil  and  have  presented  a traction 
model  based  on  the  experimental  data.  It  is  shown  that  the  viscosity  pressure 
temperature  relationship,  in  the  high  pressure  zone,  for  the  purpose  of  com- 
puting traction  is 


y ■ 1*45  X 10”^  X (U)  exp 

,-5 


-9.949  + 


5-24719  X 10~ 
T-345 


4-18935  X 10  ' -5  v 

+ 3-60  X 10  p) 

(T-345) J 
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where 


X(U) 


1 for  U,  the  rolling  speed 


< 900  in/sec 
U > 900  in/sec 


and  T and  p are,  respectively,  the  temperature  (°R)  and  pressure  (Ibf/ln^). 
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Using  the  above  equation,  traction  Is  determined  by  simultaneously  Integrating 
the  Newtonian  shear  and  the  energy  equation  through  the  film  when  the  required 
boundary  conditions  are  determined  by  the  prescribed  surface  temperature  and 
the  rolling  speed. 


It  turns  out  that  such  an  integration  has  to  be  performed  numerically  and, 

hence,  the  model  becomes  somewhat  time  consuming.  In  view  of  this  difficulty, 

equation  (V.19)  is  further  expanded  to  make  it  similar  to  (V.16),  and 

coefficients  similar  to  u ",  ot",  and  Q*  are  derived.  These  coefficients  will 

o 

now,  of  course,  be  functions  of  temperature  and  some  straightforward  algebra 
will  show  that  the  required  expressions  are 


/"h  . ,«-12  [5-24719  X 10^ 

X(U)  6-92743  x 10  exp  J 

[ (T  - 345) 

5-24719  X 10^  _ 8-37871  x 10^  ^ ^o^^ 

(T  - 345)^  (T  - 345)^ 


4-18935  X 10^?  Ibf-sec 
(T  - 345)^  j in^ 

(V.20) 


a*  * 3-6  X 10  ^ Ibf/in^ 


x(U)  is  the  same  as  that  described  above  in  equation  (V.19),  and  T will  be 
the  inlet  temperature  (°R) . 


If  the  above  values  are  used,  then  the  7808  model  becomes  identical  to  the 

* * 

polyphenyl  ether  model,  except  that  the  apparent  coefficients,  , a , and 
3*  are  temperature  dependent.  The  differences  between  such  a simplified 
model  and  the  actual  model  (reference  13)  are  shown  in  Figures  V-3  and  V-4 
for  some  extreme  cases.  It  is  clear  from  these  comparisons  that  the  above 
approximations  are  quite  reasonable. 


This  simplified  approach  also  suggests  that  the  general  analytical  approach 
can  be  used  somewhat  universally  for  any  lubricant.  It  will  be  necessary  to 
derive  an  equation  similar  to  (V.19)  and  then,  depending  on  the  actual  form 
of  the  equation,  appropriate  expansions  may  be  substituted  and  relevant  expres- 
sions for  the  apparent  coefficients,  similar  to  equations  (v.20),  may  be 
derived.  Once  these  coefficients  are  determined,  the  remaining  analysis  will 
be  identical  to  that  described  in  the  case  of  the  polyphenyl  ether  model. 


i 


Once  again  the  pressure-temperature-vlacosity  relation  of  the  form  (V.18) 
is  used  for  the  computation  of  the  film  thickness.  The  low  pressure  viscosity 
data  for  the  7808  oil  has  been  found  to  fit  the  following  expressions 
(reference  14). 


1»45  X 10~^  exp 


exp  (-3-7048  In  + 24-394)  - 0.87 


Ibf-sec 

In'* 


^ 3(i)t) 

exp  fexp  (-3-7048  In  to  + 24-394)1  exp  (-3-7048  In  To  + 24-394)  3.7048  To 
exp  [exp  (-3.7048  In  To  + 24.394)]  - 0.37 


vl 


5-227  X 


0-013332 
Tq  - 388 


(V.21) 


when  T^  Is  the  reference  temperature  (*^3) - 

Thermal  conductivity  of  the  MIL-L-7308  is  found  to  be  1-2062  x lO"*"  Ibf /sec/*"’?. 
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ire  V-3.  Comparison  of  Actual  and  Simplified  Traction  Models 
At  ■ / T«  rat' 


Figure  V-4.  Comparison  of  Actual  and  Simplified  Traction  Models 
At  High  Temperature 
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i i. 


1 ' 

ROLLER 

MOTION  RESULTS 

i . 

r “ 

The  analysis  presented  In  the  preceding  section  forms  a basis  for  a roller 

i . 

bearing  performance  simulation  con^uter  program.  This  can  be  seen  by 

7 - 

examining  some  results  for  the 

roller  motion  in  a typical  high-speed  roller 

i > 

bearing.  Some  geometrical  details  for  the  particular  bearing  are  as  follows 

\ 

• « 

Bore 

6.49606  in. 

1 

Outside  Diameter 

8.858  in. 

J . 

Roller  Diameter 

0.6299  in. 

: - 

Roller  Length 

0.6299  in. 

* * 

Central  Flat  Land 

0.20  in. 

I 

Crown  Radius 

17.0  in. 

I . 

Number  of  Rollers 

32 

- - 

Fitch  Diameter 

7.75  in. 

i . 

Diametral  Clearance 

5.8  X 10  ^in. 

Roller  Corner  Radius 

0.050  in. 

- - 

All  the  details  about  the  cage 

geometry  have  been  omitted  since  the  emphasis 

is  only  on  the  roller  motion  in  the  present  investigation. 

1 

« • 

• 

T 

The  operating  conditions  assumed  to  be  imposed  on  the  bearing  are: 

! 

4 ' 

Radial  Load 

1465  Ibf 

^ ' 

Inner  Race  Speed 

13230  rpm 

• « 

Outer  Race  Speed 

0. 

m • 

i 

Inner  Race  Temperature 

650°R 

i * 

Outer  Race  Temperature 

610°R 

T ' 

i. 

Lubricant 

MIL-L-7808 

• m 

] 

• « 

In  order  to  sttidy  the  motion  of  a roller  under  the  above  condition,  a simple 

!: 
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quasi-static  type  computation  is  first  performed  to  determine  the  relative 
position  of  the  races.  As  shown  in  Figure  Vl-1,  the  races  are  then  held 
in  this  eccentric  position  and  all  the  rollers  are  removed  from  the  bearing. 
One  roller  is  placed  just  before  the  start  of  the  load  zone  with  conditions 
identical  to  those  in  the  full  bearing.  Using  these  conditions  as  initial 
conditions,  the  differential  equations  of  roller  motion  are  then  integrated 
to  simulate  the  general  motion.  Since  such  an  integration  is  performed 
numerically,  the  equations  are  properly  dimensionalized  and  the  scales  for 
the  fundamental  quantities  are  derived  from  the  bearing  geometry  and  applied 


conditions . 

Force 

= Applied  Load  = 

1465  Ibf. 

Length 

= Roller  Radius  = 

0.3149  in. 

Time 

= y/iar/Q  - 

1.7495  X 10  ^ sec. 

where  m is  the  mass  of  the  roller,  r is  the  roller  radius,  and  Q is  the 
applied  load. 

The  above  quantities  will  also  be  used  in  dimens ionalizing  most  of  the 
results  discussed  below. 

With  reference  to  Figure  Vl-1,  initially  due  to  the  centrifugal  force,  the 
roller  is  in  contact  with  the  outer  race  only.  As  it  moves  in  the  load  zone, 
it  comes  in  contact  with  the  inner  race  and  the  load  variation  as  the  roller 
moves  through  the  load  zone  is  shown  in  Figure  VI-2.  The  roller  will 
experience  radial  accelerations  due  to  this  load  cycle  and  also  the  resulting 
tractive  forces  will  contribute  to  the  orbital  and  angular  accelerations  of 
the  roller.  The  general  patterns  of  these  accelerations  are  shown  in  Figure 
VI-3.  The  high  frequency  variations  in  the  radial  and  orbital  accelerations 
correspond  to  the  roller/race  contact  spring  resonance.  In  the  load  zone  these 
will  be  both  inner  and  outer  race  contact  springs  in  parallel  while  only  the 
outer  race  spring  is  relevant  in  the  unloaded  zone.  Hence,  the  natural 
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Figure  VI-2.  Roller/Race  Contact  Load  Variation  as  the  Roller 
Travels  Through  the  Load  Zone. 
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Figure  Vl-3.  Typical  Roller  Acceleration  Profiles  as  It 
‘ Travels  Through  the  Load  Zone. 


frequency  will  be  somewhat  higher  In  the  load  zone.  These  frequencies  are 
similar  to  the  ones  identified  by  Gupta,  et.al.  [15]  for  ball  bearings. 

In  general,  the  contact  resonance  frequency  is  somewhat  dependent  on  load. 
Hence,  the  frequency  in  the  load  zone  is  not  constant  in  Figure  IV-3. 

However,  since  the  centrifugal  force  is  fairly  constant,  the  resonance 
frequency  is  constant  in  the  unloaded  zone.  On  an  average,  the  load  and 
no  load  zone  frequencies  in  Figure  VI-3  are  approximately  24  and  17  kHz, 
respective ly. 

No  skew  type  motion  of  the  roller  is  observed  in  the  present  example,  since 
the  races  are  perfectly  aligned  and  the  system  is  free  of  any  imposed  moments 
about  the  transverse  axis  of  the  roller.  In  a real  bearing  with  misalignment 
or  some  geometrical  deviations  will  result  in  roller  skew.  Skidding  of  the 
roller  primarily  means  relative  sliding  between  the  roller  and  the  races. 
Figure  VI-4  shows  the  roller/race  slip  velocities  and  the  resulting  traction 
coefficients  for  both  the  outer  and  inner  race  contacts.  It  is  seen  that  the 
roller  starts  out  with  a zero  slip  (the  selected  initial  condition)  and  it 
experiences  a slip  cycle  as  it  travels  through  the  load  zone.  Since  the 
magnitude  of  slip  is  quite  small,  the  traction  curve  is  almost  linear  with 
slip  and,  therefore,  the  traction  coefficients  in  Figure  VI-4  appear  to  be 
almost  proportional  to  slip. 

These  general  results  primarily  summarize  the  practical  significance  of  the 
roller  motion  analysis,  presented  in  earlier  sections,  in  simulating  the 
dynamic  performance  of  a roller  bearing.  In  a complete  bearing  such  an 
analysis  will  be  performed  for  every  roller  and  also  a treatment  of 
roller/cage  and  race/cage  Interactions  will  be  Included. 
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Figure  VI-4.  Roller /Race  Slip  and  Traction  Variations. 
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SUMMARY  I 

j 

A generalized  formulation  of  the  differential  equations  of  motion  of  a roller 
In  a radially  loaded  cylindrical  roller  bearing  Is  presented.  The  motion  Is 
considered  In  two  parts,  e.g.,  motion  of  the  roller  mass  center  In  an  Inertial 
fr£une  of  reference  and  the  angular  motion  about  the  roller  mass  center  In  a 
roller  fixed  coordinate  frame.  The  formulation,  therefore,  has  the  complete 
six  degrees  of  freedom  and  It  has  the  capabilities  of  treating  roller  skew  and  ! 

other  complicated,  and  often  undesirable,  motion.  Also,  the  geometric  formu- 
lation takes  full  account  of  any  misalignment  of  the  races.  ^ 

I 

i 

Analytical  frame  work  for  the  computation  of  applied  forces  and  moments  at  the 
roller-race  Interaction  Is  described.  The  normal  contact  force  Is  primarily  | 

determined  by  locating  the  geometric  center  of  the  roller  with  respect  to  the 
Interacting  surface  of  the  roller  and,  therefore,  by  computing  the  elastic  deflec-  ] 

tlon  at  the  contact  point.  Knowing  the  deflection,  the  normal  contact  load  Is  ' 

determined  by  the  most  commonly  used  Falmgren  or  Lundberg  type  of  load  deflection 
relation  for  a line  contact.  The  tractive  forces  ^md  moments  are  computed  by  | 

first  determining  the  local  slip  velocities  In  the  contact  zone  and  then  estl- 
mating  the  traction  coefficient  from  the  given  tractlon-sllp  model.  The  total  H 

contact  zone  Is  divided  Into  several  elements  and  the  forces  and  moments  are  j 

computed  for  each  element  and  a proper  summation  Is  used  to  compute  the  total  h 

Interaction.  Thus,  a full  account  for  the  local  geometry,  such  as  partial  crowns. 

Is  taken. 

In  the  case  of  roller  skew  the  Interaction  between  the  roller  corner  and  the  | ^ 

race  flange  Is  considered.  The  general  approach  Is  very  similar  to  that  for  the  ! ; 

roller  and  race  Interaction  described  above.  The  position  of  the  roller  Is  | 

located  with  respect  to  the  race  flange  and  the  elastic  deflection  Is  determined  I 

by  the  geometric  Interaction.  VAien  no  contact  between  the  roller  and  flange  ' 

I ' ; 

exists,  any  hydrodynamic  Interaction  Is  neglected,  hence,  only  metal  to  metal 

contact  Is  considered.  The  tractive  forces  and  moments  are  determined  by  first 

computing  the  slip  velocity  and  then  using  a specified  traction  slip  relation.  1 | 

i;  ‘ 

ij  ■■  y 


VII-1 


Elastohydrodynamlc  traction  models  for  three  lubricants,  the  Shell  turbo-33 
oil,  5P4E  polyphenyl  ether,  and  the  MIL-L-7806  oil,  are  presented  In  a form 
which  Is  readily  adaptable  In  the  computerized  simulations.  All  models  are 
seml-emplrlcal  In  nature  and  the  characteristic  coef ''Iclents  are  derived  by 
curve  fitting  actual  traction  data  obtained  for  a wide  range  of  operating 
conditions.  For  bearings  which  have  limited  lubricant,  relationships  to 
modify  the  film  thickness  for  starved  conditions  are  presented  and,  hence, the 
the  Influence  of  starvation  on  the  traction  characteristics,  and  ultimately 
the  bearing  dynamic  behavior.  Is  treated  to  a first  approximation. 

For  a typical  bearing  the  equations  of  motion  are  numerically  Integrated  In 
order  to  examine  the  general  nature  of  roller  motion.  It  Is  shown  that  the 
radial  acceleration  of  the  roller  mass  center  demonstrate  the  presence  of  a 
roller/race  elastic  contact  resonance.  This  high  frequency  vibration  results 
In  variations  In  normal  contact  load;  It  also  leads  to  corresponding  variations 
In  the  tractive  forces  and,  hence,  the  orbital  accelerations.  The  general 
nature  of  roller/race  slip  variations,  and  the  resulting  traction  coefficients, 
as  the  roller  travels  through  the  load  zone  Is  also  simulated. 
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SECTION  VIII 


RECOMMENDATIONS  FOR  FUTURE  RESEARCH 


One  of  the  primary  findings  of  the  present  research  has  been  that  the 
roller  Is  subjected  to  a high  frequency  motion  resulting  from  roller/race 
resonance.  Such  a behavior  should  be  confirmed  experimentally.  It  Is 
expected  that  these  high  frequencies  will  generally  cause  "ringing"  of  the 
outer  race  and,  hence,  their  presence  can  be  detected  by  picking  up 
acceleration  signals  from  the  stationary  race. 

Roller  skew  and  bearing  misalignment  are  Important  subjects  which  can  be 
studied  by  using  the  computer  programs  developed  on  the  basis  of  the  general 
roller  motion  presented  In  this  report.  Such  a parometrlc  study  will  be  of 
a great  practical  and  design  significance. 
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